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RAREPIED-GAS HEAT TRANSFER BETWEEN PARALLEL FIATES BY A MONTE CARLO METHOD 
by Morris Perlniutter 
Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio 

ABSTRACT 

The problem of rarefied-gas heat transfer between parallel plates is treated by a Monte Carlo pro- 
cedure. Hard-sphere molecule collisions are assumed. The distributions of target molecules are assumed to- 
be locally Maxwellian, with thermal velocity and density varying across the channel. The temperature, 
density, and heat transfer are found for flow conditions ranging from Khudsen numbers of 0.05 to the free- 
molecule regime for a wall- temperature ratio of 1 to 4. The present results were compared with the results 
of another Monte Carlo solution obtained by using a different method and were found to agree with those 
results; they were also compared with other analytical solutions. The present results also agreed with the 
results obtained from a continuum solution with slip boundaries except that the continuum conductivity was 
too high for the larger Knudsen number cases. 


AUSZIXJ 

Die Frage der Warmeubertragung in verdunntem Gas zwischen parallelen Flatten wird mit Hilfe einer 
Monte -Carlo Methode behandelt. Zusammenstysse von Hart sphere nmolekulen werden angenommen. Es wird 
vorausgesetzt dass die Beschussmolekule ortlich eine Maxwellverteilung haben und dass die themische 
Geschwindigkeit und Dichte im Querschnitt veranderlich sind. Die Methode liefert die Teraperatur, Dichte 
und Warmeubertragung fiir Stromungsbedingungen, die zwischen Knudsenzahlen von 0.05 und dem freien 
Molekulfluss fur ein Verhaltnis von Wand zu Temperatur bei 1 : 4 liegen. Die vorliegenden Ergebnisse 
werden mit den Resultaten einer anderen Monte- Car lolHsung unter Verwendung einer unterschiedlichen Methode 
verglichen und stiramten mit die sen Resultaten uberein; sie warden auch noch mit anderen analytischen Los- 
ungen verglichen. Die gegenwsrtigen Resultate waren auch ira Gleichklang mit durch eine Kontinuumlosung 
mit Gleitlinien erhaltenen Ergebnisse aber die Kontinuumleitfahigkeit war zu hoch in Fallen mit grosseren 
Knudsenzahlen. 


AHHOTAUHH 

HccJie^yeTCH Tenjionepe^iawa paspemeHHbiM rasoM MeK;^y napajiejibHbiMH njiacTHnaMH cnoco6oM 
MoHTe-KapJio. IIpeflnojioraeTCH coy;^apeHne 3K6ctkhx c(|)epHwecKnx MOJieKyji. Hphhhto MecTHoe MaKc- 
BejiJiOBCKOe pacnpejnejieHHe MOJieKyji MHmeHH c nepeMennofi CKopocTbK) TenjiOBoro ;;BH3iceHHH h 
nepeMeHHOft njiOTHOCTbro nonepeK Kanajia. TeMnepaTypa, nJiOTHocTb h Ten^nonepe^awa HaiifleHH 
peacHMa oOTeKanHHH c wHCJiaMH Kny^cena HawHHan ot 0.05 AO peKHMa oGTeKaHHH co cboGo^humh 
MOJieKyjiaMH npH snaweHHHx K03$$Hn;neHTa "cTeHa-TeMnepaTypa” ot 1 ;^o 4 . PesyjibTaTH 6ujih cpas- 
H6HH c peayjibTaTaMH pemeHUH nojiywennoro ^pyrHM chocoGom MoHTe-Kapjio h oxasajiHCb c hhmh b 
corjiacHK. PesyjibTaTw TaK^Ke 6hjih cpaBHenM c ;^pyrHMH anajiHTHwecKHMH pemeHHHMH. Ohh Tanme b 
corjiacHH c pesyjibTaTaMH noJiyweHHUMH ajih cnJiomHoft cpeAti c rpanwhaMK CKOJibxeHMH; ho npoBO^H- 
MOCTb cn.noniHoii cpe^H OKasaJiacb wpesMepHO BwcoKofi cjiyuaeB c 6ojiee bhcokhmh uHCJiaMH 
Knyflcena. 


INTRODUCTION 


In rarefied-gas transport problems, the 
Navier-Stokes ' equation for momentum transfer or 
the Fourier equation for conduction are no longer 
applicable because these equations assume local 
isotropy and small gradients compared with the 
path lengths of the molecules. Therefore, the 
more fundamental Boltzmann equation must be re- 
sorted to. This equation is difficult to treat by 
the usual analytical procedures because of its 
complexity. 

Most of the analyses of rarefied-gas heat 
transfer have been for the linearized problem. As 
yet, few nonlinear results have been published. 

The present analysis of heat transfer between par- 
allel plates at different temperatures enclosing a 
rarefied gas with hard-sphere molecule collisions 
is by a Monte Carlo procedure. This method re- 
duces the complexity of the analysis thereby 
avoiding many of the simplifying assumptions gen- 
erally made in the usual analytical procedures but 
at the expense of added numerical computations. 



A linearized solution to the problem of heat 
transfer between peirallel plates enclosing a hard- 
sphere gas was analyzed by Gross and Ziering [l].^ 
They assume the difference in wall temperature was 
small compared to the magnitude of the wall temper- 
atures so that the problem can be described by the 
linearized Boltzmann equation. The distribution 
function was then approximated by half -range poly- 
nomials in velocity space and the space-dependent 
coefficient were determined by forming half -range 
moment equations. 

The nonlinear problem was treated by Lavin 
and Ifeviland [2] for hard -sphere molecule colli- 
sions by using Lee's moment method. This method 
consists of choosing a distribution of molecular 
velocities in the form of half Maxwellians that 
contain a number of functions related to the local 
density and thermal velocity. A sufficient number 
of moment equations are then solved to determine 
these functions. Haviland and Lavin have also 
treated this problem by a Monte Carlo procedure 


^Numbers in brackets denote references. 
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[3], [4], and [5]. The major difference "between 
the present analysis and that of Haviland and 
Lavin is in the methods of choosing target- 
molecule collision partners and evaluating the 
local mean properties. 


The Monte Carlo Procedure 


The Monte Carlo procedure is a model sampling 
technique. A model is created and histories of 
sample molecules are followed throu^ this model 
hy making choices at points of decision from the 
appropriate prohahility distribution. By averag- 
ing certain properties of the sample molecules at 
various positions, macroscopic quantities of in- 
terest can he obtained. 

The model is shown in Figure 1. The wall at 
X2 of 0 is considered to he at temperature 
T^ Q and the other wall at xg of D is at tem- 
perature T^ A rarefied, hard-sphere molecule 
gas is contained between the walls. 

The sample molecule history is started at 
wall 0 hy picking its velocity components from 
the appropriate distribution of velocities of the 
absorbed and reemitted molecules. The space be- 
tween the walls is divided into zones as shown in 
Figure 1. The sample molecule, after leaving the 
wall, will either x»-ss throu^ the first zone or 
will have a collision with a target molecule in 
this zone. This will depend on whether the path 
length to collision for the sample molecule is 
longer or shorter than the distance the sample 
molecule must travel to pass through the zone. If 
the sample molecule passes throu^ the first zone 
with no collision, it is started again at the be- 
ginning of the next zone with its velocity compo- 
nents unchanged. 

If there is a collision in the zone, however, 
the point of collision is found at the end of the 
path length to collision. The distribution of 
target -molecule velocities is assumed to have the 
form of a Majcwellian distribution throughout each 
zone. The distribution of target -molecule colli- 
sion partners is derived from this distribution of 
target molecules. A target molecule is picked from 
the distribution of target-molecule collision part- 
ners in the zone. A collision calculation is car- 
ried out for the sample- and target-molecule col- 
lision partners, and new velocity components are 
found for the sample molecules after collision. 

The history of the sample molecule is then con- 
tinued from the point of collision with the new 
sample-molecule velocity components. The collision 
partner is ignored after collision for reasons dis- 
cussed in [4]. 

A new path length to collision is found for 
the sample molecule, and this is compared in length 
with the distance from the point of collision to 
the zone boundary. If the path length to collision 
is greater, the sample molecule is started at the 
beginning of the next zone with its velocity compo- 
nents unchanged. In the other case, there is a 
collision in the zone and the collision calcula- 
tions are repeated as before. 

If the sample molecule strikes the wall at 
X2 = D, it is absorbed and reemitted with new ve- 
locity components picked from the appropriate dis- 
tribution of velocities based on the temperature of 
the wall. After leaving the upper wall it is fol- 
lowed as before. 

The transport and fluid characteristics of in- 
terest are density distribution, temperat\ire dis- 
tribution, and heat transfer across the channel. 


These can be obtained [6] by locating scoring posi- 
tions at various locations across the channel (see 
Fig. 1). By scoring the various characteristics of 
the sample molecules as they pass the scoring cross 
section, the transport properties and fluid charac- 
teristics of interest can be obtained at these 
scoring positions. 

In the present analysis it was assumed that 
the target molecules in each zone were in a 
Maxwellian velocity distribution based on the local 
temperature and density of each zone. From the 
Monte Carlo solution based on this assumed distri- 
bution of target molecules, the temperature and 
density at each scoring position are found. The 
Monte Carlo solution is then rerun with the new 
local temperature and density in the Maxwellian 
distribution of target molecules. The problem was 
iterated in this manner until the density and tem- 
perature distribution found from the sample mole- 
cule histories agreed with the density and temper- 
ature distribution used for the target molecules. 

In [3], [4], and [5], a Monte Carlo solution 
is used to treat the present problem. That solu- 
tion is carried out with an assumed distribution of 
target molecules in a tabular form. By scoring the 
velocity components of many sample molecules as 
they pass through each zone, the velocity distribu- 
tion of the molecules in each zone can be obtained 
in tabular form. These velocity distributions eu-e 
then used as the target-molecule distributions in 
the next iteration. This process is continued un- 
til the sample-molecule distribution found agrees 
with the target -molecule distribution assumed. 

Then, from the distributions in each zone, the 
macroscopic quantities of interest can be found by 
numerical integration of the moments of the dis- 
tribution. This was carried out for one hard- 
sphere -molecule case with a Knuds en number of 2 
and a wall-temperature ratio of 1 to 4. In the 
present method, it is not necessary to find and 
store the entire distributions of target molecules 
in each increment since the form of the target 
molecule distribution is assumed. 

In the present Monte Carlo solution, since the 
collision partners are ignored, only conservation 
of mass is satisfied exactly since molecules are 
not allowed to ’’disappear" in transit between the 
walls. Conservation of momentum and energy will 
only be satisfied if the correct target molecule 
distribution is assumed. 


ANALYSIS 


Velocity Components of Emitted Sample Molecule 


The sample molecule history begins at the sur- 
face at temperature T^ q. It is assumed that the 
molecules incident on tie surface are perfectly 
accommodated; that is, they are in a Maxwellian 
distribution at the wall, based on the wall tem- 
perature. This assumption is discussed in [?]. 

In this case, the normalized Maxwellian velocity 
distribution of the molecules moving away from the 
wall V2 > 0 is given as 





2P* 




r / \i 


expf— ^1 

U.q/J 

dVidVgdVj 


( 1 ) 

V2>0 


The q/m is the number density of molecules 
moving 4way from the wall, and C^^^q is the ther- 
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mal velocity based on the wall temperature 

[(2kAl)Tw,o]^^^- 

If we consider the positive xg direction 
normal to the wall, the distribution of velocities 
of the molecules leaving the wall per unit time 
per unit area is Vgf+ o [6]). This dis- 

tribution can be transformed to cylindrical coor- 
dinates = Vj. cos 0, Vg = Vg, ^3 = Vp sin 0 
and can be normalized by (Vg) = 
to give 


The Vo, Vj., and © picked give the direc- 
tion and velocity of the sample molecule as it 
leaves the wall. This same result was used in [4] 
and [7]. Finally, V^^ ^3 obtained 

from 

V3 = Vj. sin 0; = V^. cos 0 (7) 


f(e.V2,Vr) = ^ 




2 Vo 


,0 


itC 


v,0 


\Cw,0; 


Vr d0 dVg dVj, = ^^2 ^^r 


( 2 ) 


The distributions for f0fy2fy^ can be written 
separately as 


fe de = ^ (3a) 



The velocity components of our sample molecule 
leaving the surface must be picked from these dis- 
tributions. A convenient way of picking from a 
distribution for the high-speed computer is to 
transform the distribution to a uniform distribu- 
tion in R by setting the random number R equal 
to the cumulative distribution function. For 
example , 

R =jf° fe de = ± (4) 


Calculation of the Path Length A to Collision 

The probability that a sample molecule will 
collide in the incremental path length A to 
A + dA is given in [q] as 

= ( 8 ) 

where Ag is the mean free path to collision of 
the sample molecule moving at velocity Vg in 
that zone. 

A path length to collision for the sample 
molecule from this distribution can then be chosen 
by the same procedure described earlier, 

A = -Ag In R^ (9) 

To use this relation, the mean free path Ag in 
the zone must first be calculated. 

As shown in Appendix A (Eq. (A14)), the mean 
free path of a sample molecule moving at velocity 
Vg through a Maxwellian gas at density p with a 
thermal velocity C is 


Then, by using a high-speed computer to gen- 
erate a random number R between 0 and 1, 0 can 
be obtained from Equation (4) such that for a 
large number of samples picked in this manner, the 
distribution in ii^uation (3a) will be satisfied. 

Similarly, Vg can be picked from 


As _ 2^/Ws 


D 


Pa 

„1/R V 2^;/J 


( 10 ) 



Or, since picking R is equivalent to picking 
1 - R, 

Vr = (-0^,0 m ( 6 a) 

so that by picking Ry Vg can be obtained from 
Equation (6). ^ 

The velocity Vj. is obtained in the same 
manner and is given by 

^ (6b) 


where pg is the velocity of the sample molecule 
nondimens ionalized by the thermal velocity of the 
zone (us = Vs/C). Equation (lO) can then be used 
in Equation (9) to obtain the distance to colli- 
sion for the sample molecule in the zone. 


Calculating New Sample-Molecule Velocity 
Components After Collision 

The distribution of velocities of the tsirget 
molecules colliding with the sample molecule is 
given in Appendix B follows the derivation in [8]. 
From this distribution, the velocity components of 
the target-molecule collision partner as given by 
Equation (B12) are chosen. The new velocity com- 
ponents of the sample molecule after collision 
with the target -molecule collision partner are 
given, following the derivation in [4]. 


4 



* 1 
Vis =1 

/ 2\V2 

(Vls + Vit) + Vr(i - b ) H 

(lla) 


* 1 
^2s =1 

(72s + V2t) + 1 Vr(i - 2b2) 

(11b) 

♦ 

* T 

V3s=| 

/ 2\1/2 

(V3S 4 Vjt) + Vr(i - b ; E 

(11c) 

, 

where 



OJ 

00 

00 

K) 

1 

It 

Oi 

1 

Vlt)^+ (V2s- V2t)^ 

+ (V3s-V3t)y^^ 

(lid) 


The values 
two random 
equations 

of H and H are obtained by picking 
numbers that are used in the following 


H= 2Rjj- 1 

E=2%-1 and = 

(12) 


where b^ 

must be less than 1. If b^ is 



greater than 1, a new set of random numbers must 
he chosen to find H and H. 


Macroscopic Flow Properties 


The number of sample molecules that pass the 
scoring cross section at p in the positive Vg 
direction, p, divided by the number of times 

the sample molecule leaves wall 0,N can be re- 
lated to the mass flux passing in the positive Vg 
direction at p by 




P +, 0^^^+,0 


(17) 


where 0^^^+ 0' rosiss flux leaving wall 0, 
is equal io Combining Equa- 
tions (16) and (17) results in 



If Q is taken as 1^2^ Equation (18) becomes 



The average density in the channel is then ob- 
tained by averaging the density of all the scoring 
cross sections from p = 0 to P = P- 


The macroscopic fluid characteristics needed 
are density, temperature, and heat transfer across 
the channel. Scoring positions are located at 
various distances across the channel; as shown in 
Figure 1. The average quantity of Q transported 
across the scoring cross section p in the posi- 
tive xg direction, can be written as 


(M 








^C2>*)p 


(13) 


(20a) 

2 \P+,0 P+,0/ ^-»P+,oJ 

P=1 

As shown in [4], the continuum solution, Kn -► 0, 
is given by 


Pa , 1 

P+,0 Pf 


^ = 


1 + r^/^ + r^ 


"A ' 3[l - (l - 


(20b) 


where S^. p is the number of sample molecules 
passing across the scoring cross section p in 
the positive xg direction and Q is some quan- 
tity that each sample molecule carries. Simi- 
larly, the average quantity Q transported in the 
negative Xg direction is 



(p.<av^.)p 


(14) 


Since there is no net flow across the channel 


S+^p - (15a) 

and 

(<^4/V+d^)p = -(p./v2f.d5v)^ (15b) 


where 

The temperature can be obtained in a similar 
manner. The undirected kinetic energy per unit 
mass of a monatomic gas is equal to l/2(kT/M) 
times the degrees of freedom 

II T = i <V2> (21) 

From Equation (18), this can be written as 


fS+ s. 1 



where p/(p^ q) is given by Equation (19). As 
shown in [4], the collisionless solution, Kn -+ «, 
is given by 


Hence, 


P+,pC^2>+,P 



(16) 



(21a) 


and the continuum solution, Kn 0, is given by 
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(21b) 


Finally, the heat transer across the channel 
can be obtained from 


1 = |<’2(i +’!)> 



This can be nondimensionalized by 
I* ^)X]+,0 ” then, 

by using the p+/pa previously obtained, the fol- 
lowing is found: 


1/2 





(23) 


The heat transfer 
channel as in the 
The heat transfer 
result as derived 


is then averaged across the 
case of the density (Eq. (20b)). 
is compared to the collisionless 
in [4]; 



and to the continuum solution as given by 



(24b) 


RESULTS AND CONCHJSIONS 

The teirmerature »n<3 density distribution re- 
sults are shown in Figure 2 for a Knudsen number 
of 2 and wall-temperature ratios of 0.25 and 0.75. 
These results are compared with the one case car- 
ried out by the Monte Carlo procedure of [4] and 
are in good agreement. The nonlinear moment solu- 
tion results [2] have a significantly different 
slope than the Monte Carlo results. 

The linearized solution [l] is also compared 
with the present results. For the smaller wall- 
temperatxire difference case (IV^i/T^ q = 0.75), 

the results agree with the linearized solution. 

For the larger wall-temperature differences, how- 
ever, where T^^i/T^^q = 0.25, the temperature is 
different in magnitude from the present result, 
altho\i^ the shapes of the curve agree. In the 
linearized solution the centerline temperature is 
a parameter taken midway between the wall temper- 
atures. The present results have their centerline 
temperat\ire significantly below this value for the 
Knudsen number of 2; hence the large disagreement 
between the results. 


Also shown in Figure 2 are the slip continuum 
results. The . continuum equations are solved by 
using the fluid temperature (slip temperatiire) 
near the walls rather than the wall temperatures 
as the boundary temperature. By using the slip 
temperatures that give the best fit to both the 
temperature and density results from the present 
solution, the results shown are obtained. These 
results indicate that the slip continuum densities 
and temperature ciirves agree well with the present 
results. The present results show the magnitude 
of the temperature slip between the wall and the 
fluid near the wall as well as the Knudsen layer 
where the fluid temperature near the wall deviates 
from the continuum results toward the wall temper- 
ature. 

In Figure 3, the temperature and density re- 
sults for a range of Knudsen numbers are shown. 
These results are for a fixed value wall- 
temperature ratio of 0.25. They agree with the 
collisionless solution for the very large Knudsen 
number and approach the continuum solution for the 
values of the smaller Knudsen number. Figure 3 
also shows the temperature and density distribu- 
tion obtained from the aforementioned continuum 
slip solutions. 

In Figure 4, the heat flux nondimensionalized 
by the continuum heat flux is plotted as a func- 
tion of the inverse Knudsen number for the wall- 
temperature ratio of 0.25. The Monte Carlo re- 
sult of [4) was solved for the inverse Knudsen 
number of 0.5 and found to agree with the present 
results and with the linearized solution of [l]. 

The nonlinear solution gave results above the lin- 
ear solution, while the present results fell below 
the linear solution for small Knudsen numbers. 

Also shown in Figure 4 is the heat transfer 
from the continuum solution with the slip boundary 
conditions. The slip boundary heat transfer is 
very high for larger Knudsen numbers, so that the 
slip approximation does not apply near the colli- 
sionless region. This could be because the noncon- 
tinuum conductivity, which is the first term in 
the series expansion of the second approximation 
in the Enskog method [4], is not appropriate in 
the present noncontinuum case. The results of 
Figure 4 indicate that if a conductivity lower 
than the continuum conductivity is used in the 
continuum solution with slip, the results would 
be in good agreement with the present case. 

In conclusion, the present result agrees with 
the results obtained by the Monte Carlo method 
used in [4) for the one case that was carried out 
(Kh = 2, Tv^i/Tv 0 *= 0.25). The linear and the 
nonlinear solutions disagreed somewhat with the 
present result. This could be due partly to the 
approximation made in the other analysis and 
partly to the present assumption that the target 
molecules locally are in a Maxwellian distribution. 
The assumption of a local Maxwellian distribution 
for the target molecules would be most applicable 
near equilibrium conditions, that is, cases with 
small temperature gradients and short mean free 
paths. The present results are being extended to 
the more realistic two-sided Maxwellian as was 
used in [2]. 

The results indicate the continuum solution 
with slip boundaries, if the slip can be properly 
predicted, would give good results, except that 
for the larger Knudsen numbers, the continuum con- 
ductivity value is too hi^. 

The Monte Carlo method seems to give reason- 
able results for this problem. There are certain 
drawbacks to the method, however. The error is 
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reduced by running a large number of samples, 
which run into large amounts of computing time. 

The computing time needed for a sample molecule 
was found to increase exponentially as the Knuds en 
number approached zero. On the IBM 7094, typical 
sample run times for 10,000 samples were: Kn 5, 

5 minutes; Kn = 2, 6 minutes; Kh := 0.5, 10 min- 
utes; Kn = 0.2, 15 minutes; Kn =0.1, 25 minutes; 
and Kh = 0.05, 45 minutes. However, for problems 
of the present type, which are difficult to solve 
by other methods without over simplifying assump- 
tions, the present method proved very useful. 


NOMETJCLATURE 


c 

average thermal velocity, (2kT/M)^/^ 

D 

Channel height 


erf(x) 

error function (2/Vit) / 

•A 

' -u2 

e du 

f 

probability distribution function 


probability distribution function of 


molecules moving in positive, nega- 
tive V 2 direction 

Kh 

Knuds en number Ajj,/D = U/y/2 Sp^D 

k 

Boltzmann constant 


M 

mass of molecule 


N 

number of sample molecules leaving sur- 


face 0 in Monte Carlo run, propor- 


tional to flux of molecules leaving 
surface 0 

P 

zone or scoring-position number 

Pf 

last scoring-position number 

Q 

property of sample molecule 


<Q> 

averaged quantity ^ Cifd^ 


<Q>+ 

averaged quantity Qf^d\ 

V2>0 



T 

V 

Vr 

V* 


X]_,X2,X3 


random number between 0 and 1 
wall temperature ratio lAw 0 
mutual collision cross section 
number of sample molecules through 
scoring position in positive, nega- 
tive Xg direction 

absolute temperature 
molecular velocity \'l/2 

radial velocity; (vf + v|) 

velocity after collision 
coordinates 


r 

n 

© 

A 

\ 

^^R,A 

P 

a 


defined by Eqs. (B6) and (B7) 
defined by Eqs. (B6) and (B7) 
collision rate of sample molecule with 
target molecules 
path length to collision 
mean free path length 

dimensionless velocity, v/C 
defined by Eq. (All) 

mass density 

diameter of hard-sphere molecule 
angle between sample-molecule and 
target -molecule velocities 


Subscripts : 

A averaged 

p increment number 

Pj. last increment 

R relative velocity 

s sample molecule 

t target molecule 

w,0;w,l at wall 0,1 

0,1 evaluated next to wall 0,1 

1,2,3 coordinate directions 

positive or negative direction 


APPENDIX A 


EVALUATION OF MEAN FREE PATH 


The number of collisions per unit time d© 
of a sample molecule moving at velocity Vg 
throui^ target molecules in velocity volume space 
d^V-t for hard-sphere molecules is given in [8] as 

ae = pftVRS d^v^. (Ai) 


ae 


C| 



J 


XJr (A3) 


where Vr is the velocity of the target molecules 
relative to the sample -molecule velocity before 
collision, - Vg^^, and S is the 

mutual-collision cross section where cr is 

the diameter of the molecule. The relative veloc- 
ity Vp can then be written as 


Then transforming variable by a rotation, using 
Eulerian angles as given in [9], yields 

Vit = Aij vjt (A4a) 

where 



For a MELxwellian distribution of target mole- 
cules, the collision rate can be written as 


^2s 

Vis 

- V3S 

Vis 

Vrs 

Vs 

Vrs 

Vs 

Vis 

Vzs 

- ][3s 

V2S 

Vrs 

Vs 

Vrs 

Vs 


(A4b) 
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where the inverse transformation is 




where 




"2s 


"Is 


- V' 


Is 


2s 

Vs 


"3s 


- V,« V. 


5s "Is 


- Vo^ V. 


2s "3s 


and where 

= (Vl + and V = (v? + v|)' 

Then, the following is obtained: 


1/2 


ae 




M«3/2 


exp - 


V ^ + V + V 
'^It ^2t ^'st 




aVit dVgt dV3^. 


(A5a) 


(A5b) 


(A6a) 


(A7b) 


By transforming into spherical coordinates 
^It “ ^t ® ^ 

V2t == V^ cos q) 

Vst = V^ sin 9 sin cp 
The following is obtained: 


(A7) 


d© = 


pS 




exp- ^ sin q> dcp d9 dVt (A8a) 


where 


/_? ..9 \l/2 

''R = l^''t ''s - ^''sH 


(Aob) 


This can be written as 

= 2£ 

Mjt' 


'ie = ^5^[®^(-Mt)]4t ^XR sin cp d<p d9 d^^ (A9) 


where ^ is the nondimens ionali zed velocity 
V/C. To obtain the total collision rate for the 
target molecules over all velocities. Equa- 
tion (A9) must be integrated over cp, 9, and pit- 
This can be integrated over 0 from 0 to 2jt 
and over cp from 0 to n to give 


dB( 


" j£i 72" (AlOa) 


where 


= (lit + US - dtp 


%A = 


^*s 1^6 > ^*t 

l^*t + 4/^ > l"s 


(AlOb) 


(All) 


Integrating over |i^ from 0 to oo to 
obtain the total collision frequency for a sample 
molecule moving at velocity (ig through a 
Maxwellian gas will then produce 

>■^^[^‘(-.^(-.•4)] <*1^. 

The nondimens ional mean free path is then found, 
as discussed in [8l, by dividing the sample- 
molecule velocity by the total collision rate to 
give 

\ts ^ ^s _ 

° ® ^ i^s) 

(A13) 

By using the difinition of JOiudsen number given 
by Equation (ll) the following can be written: 

“ Va/ 


+ (erf 


(A14) 


APPENDIX B 

PICKING THE TARGET-MOLECULE COLLISION PARTNER 


The velocity distribution of target molecules 
that the sample molecule will collide with as dis- 
cussed in [l] is given by 


de 1 

exp(-i4) 

2 

sin cp dcp d9 dpi^ 

9 

«l 

[exp(-n|)] 

3/2 , 

+ 3t (erf Us) 1 

(“••iSr) 


The distribution in 9 for the target molecules 
is readily seen to be 

fe d9 = |i (B2) 

Then, the 0 can be picked from this distribution 
■by 
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8 


0 = 2jtRe (B3) 

The distribution of ii,. for the target mole- 
cules is obtained from the marginal distribution 


4 

exp 

(4 


[expfui)] 

+ (erf ^3) 1 

(“••la 


(B4) 


We can then pick from this distribution as before 


t = (B5) 

^ [exp(-^xi)] + (erf Hg) 


where, when Hg > 

--2^si^|fxp(-M|)] +i^s ^H;) 

(B6) 

and 

T] = 0 

and, when 

r(us,Mt) = lr(^^s>^*s) 

and 

11 = -a[exp(-MD] l^ + 1 + ^ 

+ [2 exp (p2)J 

To find qp for the target molecule, the dis- 
tribution of target molecules as a product of a 
marginal distribution times a conditional distri- 
bution can be written: 

l!0^ = f(^)f(,|p,) (B8) 

Then, q> can be picked from the conditional dis- 
tribution 


cos 


w„. X 

r(,K) 


cp) 


2 


Picking cp for a given found as before 

from 



V'lt = ht Cp cos 0 

sin 

T 

(BUa) 


V'2t H Cp cos <p 



(Bllb) 


V’st = Mt Cp sin 0 

sin 

T 

(Bile) 

and from Equation (A4) 






^^rs 

Vls\ 

^sj 

(Bl2a) 

V2t = -V’lt^ 


^ (Bl2b) 

^3t = 0 + V-2tl 




(B12c) 
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+ Us - 2 Mt>^s K - l^sl ^ 

(m|+ h|+ 2|^Hs) - K - ^sl ^ 

After picking cp, and 0, the velocity compo- 
nents before transformation must be found. These 
are obtained from 
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Figure 1. - Model. 
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Figure 2. - Distributions across channel for Knudsen number 2. 



E-3222 




Figure 3. - Distributions for various Knudsen numbers; wall temperature ratio, 0.25. 
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Inverse Knudsen number, HK^ 

Figure 4. - Dimensionless heat flux for wall-temperature ratio of 0. 25. 
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